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TAYLOR EXPANSIONS OF SOLUTIONS OF STOCHASTIC 
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Johann Wolfgang Goethe University, Frankfurt am Main 

The solution of a parabolic stochastic partial differential equation 
(SPDE) driven by an infinite-dimensional Brownian motion is in gen- 
eral not a semi- martingale anymore and does in general not satisfy 
an Ito formula like the solution of a finite-dimensional stochastic or- 
dinary differential equation (SODE). In particular, it is not possible 
to derive stochastic Taylor expansions as for the solution of a SODE 
using an iterated application of the Ito formula. Consequently, un- 
til recently, only low order numerical approximation results for such 
a SPDE have been available. Here, the fact that the solution of a 
SPDE driven by additive noise can be interpreted in the mild sense 
with integrals involving the exponential of the dominant linear op- 
erator in the SPDE provides an alternative approach for deriving 
stochastic Taylor expansions for the solution of such a SPDE. Es- 
sentially, the exponential factor has a mollifying effect and ensures 
that all integrals take values in the Hilbert space under consideration. 
The iteration of such integrals allows us to derive stochastic Taylor 
expansions of arbitrarily high order, which are robust in the sense 
that they also hold for other types of driving noise processes such 
as fractional Brownian motion. Combinatorial concepts of trees and 
woods provide a compact formulation of the Taylor expansions. 

1. Introduction. Taylor expansions are a fundamental and repeatedly 
used method of approximation in mathematics, in particular in numeri- 
cal analysis. Although numerical schemes for ordinary differential equations 
(ODEs) are often derived in an ad hoc manner, those based on Taylor ex- 
pansions of the solution of an ODE, the Taylor schemes, provide a class of 
schemes with known convergence orders against which other schemes can be 
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compared to determine their order. An important component of such Taylor 
schemes are the iterated total derivatives of the vector field corresponding 
higher derivatives of the solution, which are obtained via the chain rule; see 

An analogous situation holds for Ito stochastic ordinary differential equa- 
tions (SODEs), except, due to the less robust nature of stochastic calculus, 
stochastic Taylor schemes instead of classical Taylor schemes are the start- 
ing point to obtain consistent higher order numerical schemes, see [29] for 
the general theory. Another important difference is that SODEs are really 
just a symbolic representation of stochastic integral equations since their 
solutions are not differentiable, so an integral version of Taylor expansions 
based on iterated application of the stochastic chain rule, the Ito formula, 
is required. Underlying this method is the fact that the solution of a SODE 
is an Ito-process or, more generally, a semi-martingale and in particular of 
finite quadratic variation. 

This approach fails, however, if a SODE is driven by an additive stochastic 
process with infinite quadratic variation such as a fractional Brownian mo- 
tion, because the Ito formula is in general no longer valid. A new method to 
derive Taylor expansions in such cases was presented in [20, 23, 28]. It uses 
the smoothness of the coefficients, but only minimal assumptions on the 
nature of the driving stochastic process. The resulting Taylor expansions 
there are thus robust with respect to assumptions concerning the driving 
stochastic process and, in particular, remain valid for other noise processes. 

A similar situation holds for stochastic partial differential equations 
(SPDEs). In this article, we consider parabolic SPDEs with additive noise 
of the form 

(1) dU t = [AU t + F(U t )]dt + BdW t , U = u , t€[0,T\, 

in a separable Hilbert space H, where A is an unbounded linear operator, F 
is a nonlinear smooth function, B is a bounded linear operator and Wt, t>0, 
is an infinite-dimensional Wiener process. (See Section 2 for a precise de- 
scription of the equation above and the assumptions, we use.) Although the 
SPDE (1) is driven by a martingale Brownian motion, the solution process 
is with respect to a reasonable state space in general not a semi-martingale 
anymore (see [10] for a clear discussion of this problem) and except of special 
cases a general Ito formula does not exist for its solution (see, e.g., [10, 37]). 
Hence, stochastic Taylor expansions for the solution of the SPDE (1) cannot 
be derived as in [29] for the solutions of finite-dimensional SODEs. Conse- 
quently, until recently, only low order numerical approximation results for 
such SPDEs have been available (except for SPDEs with spatially smooth 
noise; see, e.g., [14]). For example, the stochastic convolution of the semi- 
group generated by the Laplacian with Dirichlet boundary conditions on 
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the one-dimensional domain (0,1) and a cylindrical /-Wiener process on 
H = L 2 ((0, 1),R) has sample paths which are only — e)-H61der continu- 
ous (see Section 5.4) and previously considered approximations such as the 
linear implicit Euler scheme or the linear implicit Crank-Nicolson scheme 
are not of higher temporal order. The reason is that the infinite-dimensional 
noise process has only minimal spatial regularity and the convolution of the 
semigroup and the noise is only as smooth in time as in space. This compa- 
rable regularity in time and space is a fundamental property of the dynamics 
of semigroups; see [41] or also [9] , for example. To overcome these problems, 
we thus need to derive robust Taylor expansions for a SPDE of the form (1) 
driven by an infinite-dimensional Brownian motion. 

An idea used in [24] to derive what was called the exponential Euler 
scheme for the SPDE (1), that has a better convergence rate than hitherto 
analyzed schemes, can be exploited here. It is based on the fact that the 
SPDE (1) can be understood in the mild sense, that is as an integral equation 
of the form 

(2) U t = e At u + f e A{t ~ s) F{U s )ds+ f >('~ s ) B dW s 

Jo Jo 

for all t G [0,T] rather than as an integral equation obtained by directly 
integrating the terms of the SPDE (1). (This mild integral equation form 
of the SPDE is considered in some detail in the monograph [6], (7.1) and 
(7.3.4), and in the monograph [37], (F.0.2).) The crucial point here is that 
all integrals in the mild integral equation (2) contain the exponential factor 
e A(t-s) Q £ p era tor A, which acts in a sense as a mollifier and ensures that 
iterated versions of the terms remain in the Hilbert space H. The main idea 
of the Taylor expansions presented in this article is to use a classical Taylor 
expansion for the nonlinearity F in the mild integral equation above and then 
to replace the higher order terms recursively by Taylor expansions of lower 
orders (see Section 3). Hence, this method avoids the need of an ltd formula 
but nevertheless yields stochastic Taylor expansions of arbitrarily high order 
for the solution of the SPDE (1) (see Section 5.1 for details). Moreover, these 
Taylor expansions are robust with respect to the type of noise used and can 
easily be modified to other types of noise such as fractional Brownian motion. 

The paper is organized as follows. In the next section, we describe pre- 
cisely the SPDE that we are considering and state the assumptions that we 
require on its terms and coefficients and on the initial value. Then, in the 
third section, we sketch the idea and the notation for deriving simple Taylor 
expansions, which we develop in section four in some detail using combi- 
natorial objects, specifically stochastic trees and woods, to derive Taylor 
expansions of an arbitrarily high order. We also provide an estimate for 
the remainder terms of the Taylor expansions there. (Proofs are postponed 
to the final section.) These results are illustrated with some representative 
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examples in the fifth section. Numerical schemes based on these Taylor ex- 
pansions are presented in the sixth section. 

2. Assumptions. Fix T > and let (Q,J-,¥) be a probability space with 
a normal filtration Tt, t£ [0,T]; see, for example, [6] for details. In addition, 
let (H, (■,■)) be a separable M-Hilbert space with its norm denoted by | • | 
and consider the SPDE (1) in the mild integral equation form (2) on H, 
where Wt, t G [0, T], is a cylindrical Q- Wiener process with Q = I on another 
separable M-Hilbert space (U, (•,•)) (space-time white noise), B : U — > H is a 
bounded linear operator and the objects A, F and uo are specified through 
the following assumptions. 

Assumption 1 (Linear operator A). Let I be a countable or finite set, 
let (Aj)jez CR be a family of real numbers with infjgiAj > — oo and let 
( e i)i£i C H be an orthonormal basis of H. Assume that the linear operator 
A : D(A) C H -> H is given by 

for all v £ D(A) with D(A) = {v e H\ J2 ieI | | 2 1 <e^ , u> | 2 < oo}. 

Assumption 2 (Drift term F). The nonlinearity F : H ^ H is smooth 
and regular in the sense that F is infinitely often Frechet differentiable and 
its derivatives satisfy s\xp v&H \F^\v)\ < oo for all z G N := {1, 2, . . .}. 

Fix k > with sup igI («; + Aj) > and let D((k - A) r ), r£l, denote the 
domains of powers of the operator k — A : D(k — A) = D(A) C H — >■ H , see, 
for example, [41]. 

Assumption 3 (Stochastic convolution). There exist two real numbers 
7£ (0,1), 5 E (0, ^] and a constant C > such that 

|(k - A)T e As 5|| s < oo, f \e As B\ 2 m ds < Ct 25 

Jo 

holds for all t G [0,1], where | • |hs denotes the Hilbert-Schmidt norm for 
Hilbert-Schmidt operators from U to H . 

Assumption 4 (Initial value uo). The Fq/B{D({k — A) 7 ))-measurable 
mapping uq : — > D((k — A) 1 ) satisfies E|(k — A^uq^ < oo for every p £ 
[1, oo), where 7 G (0, 1) is given in Assumption 3. 
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Similar assumptions are used in the literature on the approximation of 
this kind of SPDEs (see, e.g., Assumption H1-H3 in [15] or see also [24, 32- 
34]). This setup also includes trace class noise (see Section 5.5) and finite- 
dimensional SODEs with additive noise (see Section 5.2). 

Henceforth, we fix to £ [0,T) and denote by V the set of all adapted 
stochastic processes 

X:Sl->C([to,T\,H) with sup \X t \ LP < oo Vp > 1, 

t <t<T 

and with continuous sample paths, where \Z\lp '■= (E|Z| p ) 1//p is the L p -norm 
of a random variable Z H. Under Assumptions 1 and 3, it is well known 
that the stochastic convolution 

f e A ^BdW s , t€[to,T\, 
Jo 

has an (up to indistinguishability) unique version with continuous sample 
paths (see Lemma 5 in Section 7.3). From now on, we fix such a version of the 
stochastic convolution. Hence, under Assumptions 1-4 it is well known that 
there is a pathwise unique adapted stochastic process U : ft — > C([0,T], H) 
with continuous sample paths, which satisfies (2) (see Theorems 7.4 and 7.6 
in [6]). Even more, this process satisfies 

(3) sup Ay U t \ L P <oo 

0<t<T 

for all p > 1. 

3. Taylor expansions. In this section, we present the notation and the 
basic idea behind the derivation of Taylor expansions. We write 

AUf.= U t -U t0 , At:=t-t 

for t € [to,T] C [0,T], thus AU denotes the stochastic process AU t , t G [t ,T], 
in V . First, we introduce some integral operators and an expression relating 
them and then we show how they can be used to derive some simple Taylor 
expansions. 

3.1. Integral operators. Let j 6 {0,1,2,1*}, where the indices {0,1,2} 
will label expressions containing only a constant value or no value of the 
SPDE solution, while 1* will label a certain integral with time dependent 
values of the SPDE solution in the integrand. Specifically, we define the 
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stochastic processes I® G V by 

' e AAt -m , 



e A ^F(U t() )ds, 
e A (ts) B dW s , 



to 



to 



,A(t-s) 



F(U 8 ) ds, 



3 = 0, 
3 = h 

J = 2, 
3 = 1*, 



for each t G [to,3"]. Note that the stochastic process jJe^-'JBdW, for 
tG [to, 31 given by 



z Mt-s) Bdw 



to 



BdW* 



o A(t-t ) 



to 



e A(to-s) Bdm 



for every t G [to ,3"] is indeed in P. Given iGN and 7 G {1, 1*}, we then define 
the i-multilinear symmetric mapping 7j : V 1 — > V by 

7j[ 5l ,... j5i ](t) :=i [* e A <. t -'W { >(U t0 )( 9l (s),...,g i (s))d8, 



to 



when j = 1 and by 
3jbi,...,ffi](t) 



to 



A(t-)( / FW(C/ t0 +rAC/ s )( 5l ( s ),..., 5i ( S ))- 



.)(*-!) 



cir ) (is, 



when j = 1* for all t G [to, 3] and all g\, . . . ,gi G V . One immediately checks 
that the stochastic processes 7? £P, j G {0,1,2,1*}, and the mappings 

I) : V i -> 7>, i G N, j G {1, 1*}, are well defined. 
The solution process £7 of (2) obviously satisfies 

(4) AU t = {e AAt - I)U t0 + f e A{t - s) F{U s )ds+ I e A ^BdW s 

J to J to 

or, in terms of the above integral operators, 

A[/ t = 7 (t) + 7 1 ,(t) + 7 2 °(t) 
for every t G [to, 3], which we can write symbolically in the space V as 

(5) A[/ = 7 ° + 7 1 °, +7 2 °. 

The stochastic processes 7j[<7i, . . . for g\, . . . ,gi G V, j G {0, 1, 2} and i G 
{0, 1,2,.. .} only depend on the solution at time t = to- These terms are there- 
fore useful approximations for the solution process Ut, t G [to ) 3 1 ]. However, 
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the stochastic processes I\* [g\, . . . ,gi\ for g\, ... ,gi G V and i E {0, 1,2, . . .} 
depend on the solution path Ut with £ G [toj^l- Therefore, we need a further 
expansion for these processes. For this, we will use the important formula 

7°, =7? + 7^ [AC/] 

(6) =i° + il4i° ] + i 1 AiU + il4i 2], 

which is an immediate consequence of integration by parts and (5), and the 
iterated formula 

l\, [ 91 , . . . , 9i ] = Ilfa, . . . , 9i ] + [AC/, 91 , . . . , 9i ] 

(7) =P 1 [g 1 ,...,g i }+4 +1) [I°,g 1 ,...,g i ] 

+ 4 +1) [I° 1 *,g 1 ,...,g i ]+4 +1) [llg 1 ,...,g i } 

for every g\ , . . . , g^ GV and every i£N (see Lemma 1 for a proof of the 
equations above). 

3.2. Derivation of simple Taylor expansions. To derive a further expan- 
sion of (5), we insert formula (6) to the stochastic process 1^* j that is, 

into (5) to obtain 

AU = 7° + (J? + /J. [7°] + /J. [7°, ] + ^ [7°] ) + 7° , 
which can also be written as 

(8) AU = 7° + 7° + 7° + /J. [7 °] + I\, [7°,] + 1{. [7 2 ]. 

If we can show that the double integral terms 7|„[7q], l\» [I®*] and 7j,[7 2 ] 
are sufficiently small (indeed, this will be done in the next section), then we 
obtain the approximation 

(9) Af7«ig + I? + l£, 

or, using the definition of the stochastic processes Iq, if and 7 2 , 

AU t « (e AAt - I)U t0 + ! e A ^F(U t0 )ds+ fe A ^BdW s 

J to J to 

for t G [io^T 1 ]- Hence, 

(10) C/ t « e AAt C/ to + Qf da) F(C/ to ) + jf * e^^T? dW s 

for i G [ioi^l is an approximation for the solution of SPDE (1). Since the 
right-hand side of (10) depends on the solution only at time to, it is the first 
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nontrivial Taylor expansion of the solution of the SPDE (1). The remainder 
terms l\* [Iq], l\* [if*] and l\* [Iq] of this approximation can be estimated by 

[I° ](t) + II. [J?.](t) + 1}. [I° 2 ](t)\ L , < C(At f^^ 

for every t G [to, T] with a constant C > and where 7 G (0, 1) and 5 G (0, h] 
are given in Assumption 3 (see Theorem 1 in the next section for details). 

We write Yt = 0((At) r ) with r > for a stochastic process Y £ V, if 
I^Il 2 — C(At) r holds for all i G [to,T] with a constant C > 0. Therefore, we 
have 

tft- (V A ^ + (^j\ As ds^F{U t0 ) + j\ A ^Bdw)j 
= 0((At) ( - 1+min ^^) 



or 



(11) 



+ (^ A * e As ds^ F(U t0 ) + jf e A < ( - s 'B^ ; 



+ 0((At) 



(l+min(7,5)) n 



The approximation (10) thus has order 1 +min(7,5) in the above strong 
sense. It plays an analogous role to the simplest strong Taylor expansion in 
[29] on which the Euler-Maruyama scheme for finite-dimensional SODEs is 
based and was in fact used in [24] to derive the exponential Euler scheme 
for the SPDE (1). Note that the Euler-Maruyama scheme in [29] approx- 
imates the solution of an SODE with additive noise locally with order |. 
Here, in the case of trace class noise, we will have 7 = 5 — 6, 5 = \ (see 
Section 5.5), and therefore the exponential Euler scheme for the SPDE (1) 
in [24] also approximates the solution locally with order I — e (see Section 
5.5.2), while other schemes in use, in particular the linear implicit Euler 
scheme or the Crank-Nicolson scheme, approximate the solution with order 
\ instead of order | as in the finite-dimensional case. Therefore, the Taylor 
approximation (11) attains the classical order of the Euler approximation for 
finite-dimensional SODEs and in general the Taylor expansion introduced 
above lead to numerical schemes for SPDEs, which converge with a higher 
order than other schemes in use (see Section 6). 



3.3. Higher order Taylor expansions. Further expansions of the remain- 
der terms in a Taylor expansion give a Taylor expansion of higher order. To 
illustrate this, we will expand the terms l\* [Iq] and l\* [I®] in (8). From (7), 
we have 



i\* [i°o] = n [i°o] + ii [ig, ig] + 4 [# , ig] + 4 [4 A] 
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and 

iui° 2 ] = n[i° 2 ] + n [j °, j 2 °] + ii [/?. , j 2 °] + /?, [ii 4\ , 

which we insert into (8) to obtain 

AU = 1° + If + 1 2 ° + A 1 [J °] + 1} [/ 2 °] + i?, 
where the remainder term R is given by 

R = Ih [i?.] + 1\, [i °, i °] + I?. [4 , i Q °] + II [i 2 °, io°] 

+ / 1 2 ,[/ M 2 ] + / 1 2 ,[i 1 *,i 2 °] + i 1 2 .[i 2 ,i 2 ]. 

From Theorem 1 in the next section, we will see R = 0((Ai)( 1+2mm ( 7,<5 ))). 
Thus, we have 

AU = / ° + 1? + J 2 ° + 1, 1 [7 °] + I 2 a + 0((At)( 1+2min ^ 5 ))), 
which can also be written as 

U t = e AAt U t0 + ^ At e As ds^jF(U t0 ) + j\ A ^BdW s 
+ f t e A ^F'(U to )(e AAs -I)U to d S 

+ t e A ^F'{U t0 ) f\ A ^BdW r ds + 0{(At^ 1+2m[n ^y). 

J to J to 

This approximation is of order 1 + 2min(7,<5). 

By iterating this idea, we can construct further Taylor expansions. In par- 
ticular, we will show in the next section how a Taylor expansion of arbitrarily 
high order can be achieved. 

4. Systematic derivation of Taylor expansions of arbitrarily high order. 

The basic mechanism for deriving a Taylor expansion for the SPDE (1) was 
explained in the previous section. We illustrate now how Taylor expansions 
of arbitrarily high order can be derived and will also estimate their remainder 
terms. For this, we will identify the terms occurring in a Taylor expansions 
by combinatorial objects, that is, trees. It is a standard tool in numerical 
analysis to describe higher order terms in a Taylor expansion via rooted 
trees (see, e.g., [2] for ODEs and [1, 38-40] for SODEs). In particular, we 
introduce a class of trees which is appropriate for our situation and show 
how the trees relate to the desired Taylor expansions. 
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4.1. Stochastic trees and woods. We begin with the definition of the trees 
that we need, adapting the standard notation of the trees used in the Taylor 
expansion of SODEs (see, e.g., Definition 2.3.1 in [38] as well as [1, 39, 40]). 

Let N E N be a natural number and let 

t':{2,..., N} -> {1,... ,N-1}, t":{l,.. .,#}-> {0,1,2,1*} 

be two mappings with the property that t' (j) < j for all j E {2, . . . , N}. The 
pair of mappings t = (t',t") is a S-tree (stochastic tree) of length N = Z(t) 
nodes. 

Every S-tree can be represented as a graph, whose nodes are given by the 
set nd(t) := {1, . . . , N} and whose arcs are described by the mapping t' in 
the sense that there is an edge from j to t'(j) for every node j E {2, . . . , N}. 
In view of a rooted tree, r' also codifies the parent and child pairings and 
is therefore often referred as son-farther mapping (see, e.g., Definition 2.1.5 
in [38]). The mapping t" is an additional labeling of the nodes with t"(j) E 
{0, 1, 2, 1*} indicating the type of node j for every j E nd(t). The left picture 
in Figure 1 corresponds to the tree ti = (t'^t") with nd(ti) = {1,2,3,4} 
given by 

ti(4) = l, ti(3) = 2, ti(2) = l 

and 

t'/(l) = l, t'/(2) = r, t'/(3) = 2, t?(4) = 0. 

The root is always presented as the lowest node. The number on the left 
of a node in Figure 1 is the number of the node of the corresponding tree. 
The type of the nodes in Figure 1 depends on the additional labeling of 
the nodes given by t'{. More precisely, we represent a node j E nd(ti) by 
® if t'l(j) = 0, by • if t'{(j) = 1, by O if t'{ij) = 2, and finally by ■ if 




Fig. 1. Two examples of stochastic trees. 
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i (g) i ■ i o 

Fig. 2. The stochastic wood wo in SW. 

t"(j) = 1*. The right picture in Figure 1 corresponds to the tree t 2 = (t 2 ,t 2 ) 
with nd(t2) = {1, . . . , 7} given by 

t 2 (7)=4, t 2 (6)=4, t' 2 (5) = l, 
t' 2 (4) = l, t' 2 (3) = l, t 2 (2) = l 

and 

t 2 '(l)=0, t 2 '(2) = 0, t 2 '(3) = 2, 4(4) = 1, 
t 2 '(5) = l*, t 2 '(6) = l, t 2 '(7) = 0. 

We denote the set of all stochastic trees by ST and will also consider a tuple 
of trees, that is, a wood. The set of S-woods (stochastic woods) is defined 
by 



SW:= (J(ST) n . 



n=l 

Of course, we have the embedding ST C SW. A simple example of an S- 
wood which will be required later is wo = (ti,t 2 ,t3) with ti, t 2 and t3 given 
by Z(ti) = l(t 2 ) = l(t 3 ) = 1 and t'/(l) = 0, t 2 '(l) = 1*, t 3 '(l) = 2. This is shown 
in Figure 2 where the left tree corresponds to ti, the middle one to t 2 and 
the right tree corresponds to t3. 

4.2. Construction of stochastic trees and woods. We define an operator 
on the set SW, that will enable us to construct an appropriate stochastic 
wood step by step. Let w = (ti, . . . ,t n ) with n G N be a S-wood with tj = 
(t.,t") G ST for every i £ {1, 2, . . . ,n}. Moreover, let % £ {1, ... ,n} and j £ 
{1, . . . , /(tj)} be given and suppose that t'-(j) = 1*, in which case we call the 
pair an active node of w. We denote the set of all active nodes of w 
by acn(w). 

Now, we introduce the trees t n+ i = (t^ +1 ,t^ +1 ), t n+2 = (t^ +2 ,t^ +2 ) and 
t n+3 = K+sX+s) in ST by nd(t n+m ) = {1, . . . , l(U), i(t<) + 1} and 

^n+m(k) = ^i(k), k = 2,...,l(ti), 

*n+m(^) = *i (^)) k = 1, . . . ,l(ti), 



t'n+miKti) + 1) = 3, < +m (l(ti) + 1) = | 



1*, m = 2, 

(m — 1), else, 
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2 ® 2 ■ 2 Q 



1® 1 • 1 O 1 ■ IB !■ 

Fig. 3. The stochastic wood wi in SW. 

for m = 1, 2, 3. Finally, we consider the S-tree = (t', t") given by t' = t£, 
but with i"(k)=t"(k) for fc ^ j and t"(j) = 1. Then, we define 

En ,iw = £V; ,-)(ti, . . . ,t n ) := (ti, . . . ,tj_i,t^ J \tj + i, . . . ,t n+3 ) 



(iJ)W = -0(^)1 

and consider the set of all woods that can be constructed by iteratively 
applying the -E'(ij) operations, that is, we define 



SW' := {w } U { w G SW 



3n G N,ii, ji, . . .,i n ,j n £N: 

VZ = l,...,n:(^ 5 i/) G acn^^^) ■ • ■ 

w = E (i n ,j n ) ■ ■ ■ E (h,ji) w 



for the wo introduced above. To illustrate these definitions, we present some 
examples using the initial stochastic wood wo given in Figure 2. We present 
these examples here in a brief way, and, later in Section 5.1, we describe 
more detailed the main advantages of the particular examples considered 
here. First, the active nodes of wo are acn(wo) = {(2,1)}, since the first 
node in the second tree in wo is the only node of type 1*. Hence, £p,i) w o 
is well defined and the resulting stochastic wood wi = -B(2,i) w 0) which has 
six trees, is presented in Figure 3. Writing wi = (ti, . . . ,te), the left tree in 
Figure 3 corresponds to ti, the second tree in Figure 3 corresponds to t2, 
and so on. Moreover, we have 



(12) 



acn( Wl ) = {(4,l),(5,l),(5,2),(6,l)} 



for the active nodes of wi, so W2 = ^(4,1) Wi is also well defined. It is 
presented in Figure 4. In Figure 5, we present the stochastic wood W3 = 
-^(6.1) w 2> which is well defined since 



1 (x) 1 



2 2 (5? 



2 ® 2 



1 O 1 • 1 ■ 1 ■ 1 ■ 1 

Fig. 4. The stochastic wood W2 in SW 
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2 (X) 2 



1 ® 1 




2 (X) 



1 ■ 1 ■ 1 

Fig. 5. The stochastic wood W3 in SW. 

(13) acn(w 2 ) = {(5, 1), (5, 2), (6, 1), (7, 1), (8, 1), (8, 3), (9, 1)}. 

For the S-wood W3, we have 



(14) 



acn(\V3) 



(5,1), (5, 2), (7,1), (8,1), (8, 3), (9,1), 
(10,1), (11,1), (11,3), (12,1) 



Since (7, 1) G acn(w3), the stochastic wood W4 = Enuxv^ is well defined and 
presented in Figure 6. For the active nodes, we obtain 



(15) acn(w4) 



(5,1), (5,2), (8,1), (8,3), (9,1), (10,1), (11,1), 
(11,3), (12,1), (13,1), (14,1), (14,4), (15,1) 



Finally, we present the stochastic wood W5 = Er 12, 1)^(10, i)-^(9,i) w 4 with 

(5, 1), (5, 2), (8, 1), (8, 3), (11, 1), (11,3), (13, 1), (14, 1), ' 
(16) acn(w 5 ) = { (14, 4), (15, 1), (16, 1), (17, 1), (17, 4), (18, 1), (19, 1), 
(20, 1), (20,4), (21, 1), (22, 1), (23, 1), (23,4), (24, 1) 




Fig. 6. The stochastic wood W4 in SW. 
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Fig. 7. The stochastic wood W5 in SW. 



in Figure 7. By definition, the S-woods wo,wi, . . . , W5 are in SW', but the 
stochastic wood given in Figure 1 is not in SW'. 

4.3. Subtrees. Let t = (t',t") be a given S-tree with Z(t) > 2. For two 
nodes k, I G nd(t) with k < I, we say that Z is a grandchild of fc if there exists a 
sequence k\ = k < < • • • < k n = I of nodes with n G N such that t'(fc„+i) = 
k v for every u € {1, . . . , n — 1}. Suppose now that j\ , . . . , j n G nd(t) with n G N 
and ji < • • • < in are the nodes of t such that t'(jj) = 1 for every i = 1, . . . , n. 
Moreover, for a given i G {1, . . . , n} suppose that j^i, • • • , G nd(t) 
with = ji t i < ji t 2 < ■ • ■ < jifo < l(t) and G N are the grandchildren of ji. 
Then, we define a tree tj = (t^,t") G ST with /(tj) := li and 

ii,t<(fc) = t'0'i,fc), t"(k) = t"(j i;k ) 

for all k G {2, . . . , ZJ and t^'(l) = t"(jj). We call the trees t x , . . . ,t n G ST 
defined in this way the subtrees of t. For example, the subtrees of the right 
tree in Figure 1 are presented in Figure 8. 
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Fig. 8. Subtrees of the right tree in Figure 1. 



4.4. Order of a tree. Later stochastic woods in SW' will represent Tay- 
lor expansions and Taylor approximations of the solution process U of the 
SPDE (1). Additionally, we will estimate the approximation orders of these 
Taylor approximations. To this end, we introduce the order of a stochastic 
tree and of a stochastic wood, which is motivated by Lemma 4 below. More 
precisely, let ord : ST — > [0, oo) be given by 

ord(t) := Z(t) + ( 7 - l)|{j G nd(t)|t"(j) = 0}| 

+ (J-l)|{jGnd(t)|t"(j)=2}| 

for every S-tree t = (t',t") G ST. For example, the order of the left tree in 
Figure 1 is 2 + 7 + 5 and the order of the right tree in Figure 1 is 3 + 37 + 
5 (since the right tree has three nodes of type 0, three nodes of type 1, 
respectively, 1*, and one node of type 2). 

In addition, we say that a tree t = (t',t") in ST is active if there is a 
j G nd(t) such that t"(j) = 1*. In that sense a S-tree is active if it has an 
active node. Moreover, we define the order of an S-wood w = (ti, . . . ,t n ) G 
SW with n G N as 

ord(w) := min{ord(tj), 1 < i < n|tj is active}. 

To illustrate this definition, we calculate the order of some stochastic woods. 
First of all, the stochastic wood in Figure 2 has order 1, since only the middle 
tree in Figure 2 is active. More precisely, the node (2, 1) of the S-wood wo 
is an active node and therefore the second tree is active. The second tree in 
Figure 2 has order 1 (since it only consists of one node of type 1*). Hence, the 
S-wood wo has order 1. Since the last three trees are active in the stochastic 
wood wi in Figure 3 [see (12) for the active nodes of Wj], we obtain that 
the stochastic wood in Figure 3 has order 1 + min(7, 5). The last three trees 
in the S-wood wi have order 1 + 7, 2 and 1 + 5, respectively. As a third 
example, we consider the S-wood W2 in Figure 4. The active nodes of W2 
are presented in (13). Hence, the last five S-trees are active. They have the 
orders 2, 1 + 5, 1 + 27, 2 + 7 and 1 + 7 + 5. The minimum of the five real 
numbers is 1 + min(27, 5). Therefore, the order of the S-wood W2 in Figure 4 
is 1 + min(27, 5). A similar calculation shows that the order of the stochastic 
wood W3 in Figure 5 is 1 + 2 min(7, 5) and that the order of the stochastic 
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wood W4 in Figure 6 is 1 + min(37,7 + 5,25). Finally, we obtain that the 
stochastic wood W5 in Figure 7 is of order 1 + 3min(7, 5, |). 

4.5. Trees and stochastic processes. To identify each tree in ST with 
a predictable stochastic process in V, we define two functions : ST — > V 
and ^ : ST — > V , recursively. For a given S-tree t = (t',t") G ST, we define 
$(t) := J t °„ (1) when t"(l) G {0,2} or Z(t) = 1 and, when Z(t) > 2 and t"(l) G 

{1, 1*}, we define 

*(t) :=I^ {1) mti),---,Htn)}, 

where ti, . . . ,t n G ST with n G N are the subtrees of t. In addition, for an 
arbitrary t G ST, we define ^(t) := if t is an active tree and ^(t) = <&(t) 
otherwise. Finally, for a S-wood w = (ti, . . . , t n ) with n G N we define <3?(w) 
and ^>(w) by 

$(w) = $(ti) + • • • + $(t n ), (w) = *(ti) + • • • + *(t„). 

As an example, we have 

(17) $(w )=ig + I?.+i§ and tf(w ) = ig + i£ 

for the elementary stochastic wood wo (see Figure 2). Hence, we obtain 

(18) $(w ) = AC7 

from (5) and (17). Since (2, 1) is an active node of wo, we obtain 

(19) *( W1 ) = I° + If + I° 2 + I{* [I °] + I\. [IU + II [/ 2 °] 
and 

(20) M/(w 1 )=/ ° + / 1 ° + / 2 

for the S-wood wi = £ ; ( 2 ,i) w o presented in Figure 3. Moreover, in view of 
(6) and (7), we have 

(21) $(w) = $(E {id) w) 

for every active node (i,j) G acn(w) and every stochastic wood w G SW'. 
Hence, we obtain 

$(wi) = $(%l)W ) = $(w ) 

due to the equation above and the definition of wi. Hence, we obtain 
<I>(wi) = AU, which can also be seen from (19), since the right-hand side 
of (19) is nothing other than (8). We also note that the right-hand side of 
(20) is just the exponential Euler approximation in (9), so we obtain 

AZ7 = $(wi) ~ *(wi). 

With the above notation and definitions we are now able to present the main 
result of this article, which is a representation formula for the solution of 
the SPDE (1) via Taylor expansions and an estimate of the remainder terms 
occurring in the Taylor expansions. 



TAYLOR EXPANSIONS FOR SPDE'S 



17 



Theorem 1. Let Assumptions 1~4 be fulfilled and let w G SW' be an 
arbitrary stochastic wood. Then, for each p G [1, oo), there is a constant C p > 
such that 



holds for every t G [io>^1; where Ut, t G [0,T], is the solution of the SPDE 
(1). Here the constant C p > is independent of t and to but depends on p 
as well as w, T and the coefficients of the SPDE (1). 

The representation of the solution here is a direct consequence of (18) 
and (21). The proof for the estimate in (22) will be given in Section 7. Here, 
$(w) = AU is the increment of the solution of the SPDE (1), while ^(w) is 
the Taylor approximation of the increment of the solution and <£(w) — 'I'(w) 
is its remainder for every arbitrary w G SW'. Since there are woods in SW' 
with arbitrarily high orders, Taylor expansions of arbitrarily high order can 
be constructed by successively applying the operator Eu^\ to the initial S- 
wood wo- Finally, the approximation result of Theorem 1 can also be written 
as 



for every stochastic wood w G SW'. Here, we also remark that Assumptions 
1-4 can be weakened. In particular, instead of Assumption 3, one can assume 
that the nonlinearity F : V — > V is only i-times Frechet differentiable with i G 
N sufficiently high and that the derivatives of F satisfy only local estimates, 
where V C H is a continuously embedded Banach space. Nevertheless, it 
is usual to present Taylor expansions for stochastic differential equations 
under such restrictive assumptions here (see [29]) and then after considering 
a particular numerical scheme one reduces these assumptions by pathwise 
localization techniques (see [11, 26] for SDEs and [22] for SPDEs). 

5. Examples. We present some examples here to illustrate the Taylor 
expansions introduced above. 

5.1. Abstract examples of the Taylor expansions. We begin with some 
abstract examples of the Taylor expansions. 

5.1.1. Taylor expansion of order 1. The first Taylor expansion of the 
solution is given by the initial stochastic wood wo (see Figure 2), that is, we 
have $(wo) = AU approximated by ^(wo) with order ord(wo). Precisely, 
we have 



(22) 



U t = U t0 + <f>(w)(t), 
(E[\U t - U t0 - *(w)(t)H) 1/p <C p -(t- t ) ord(w) 



U t = U t0 +*( W ) + O((Atr d W) 
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and 

*(w )(t) = (e AAt - I)U t0 + f e A ^F(U s )ds + f e A ^BdW s 

Jto Jto 

for every t G [io>2l due to (17). Since ord(wo) = 1 (see Section 4.4), we 
finally obtain 

U t = e AAt U t0 + I e A(t ~^BdW s + O(At) 
Jt 

for the Taylor expansion corresponding to the S-wood wo. 

5.1.2. Taylor expansion of order l+min(7, S). In order to derive a higher 
order Taylor expansion, we expand the stochastic wood wo- To this end, 
we consider the Taylor expansion given by the S-wood wi = i)Wq (see 
Figure 3). Here, $(wi) and $(wi) are presented in (19) and (20). Since 
ord(wi) = 1 + min(7, 5) (see Section 4.4), we obtain 

U t = e AAt U t0 + Qf * e As ds^j F(U t0 ) + J* e A ^BdW s + 0((At) ( 1+mi ^ 6 » ) 

for the Taylor expansion corresponding to the S-wood wi. This example 
corresponds to the exponential Euler scheme introduced in [24], which was 
already discussed in Section 3.2 [see (11)]. 

5.1.3. Taylor expansion of order 1 + min(27, 5). For a higher order Tay- 
lor expansion, we now have several possibilities to further expand the stochas- 
tic wood wi. For instance, we could consider the Taylor expansion given by 
the stochastic wood ^(5,1)^(5,2) w i [ see (12) for the active nodes of wj. 
Since our main goal is always to obtain higher order approximations with 
the least possible terms and since the fifth tree of wi [given by the nodes 
(5,1) and (5,2)] is of order 2 (see Section 4.4 for details), we concentrate 
on expanding the lower order trees of Wi . Finally, since oftentimes 7 < 5 in 
examples (see Section 5.4 and also Section 5.5), we consider the stochastic 
wood W2 = i?(4 i)Wi (see Figure 4). It is of order 1 + min(27, S) (see Section 
4.4) and the corresponding Taylor approximation W(w2) of $(w2) = AC/ is 
given by *(w 2 ) = i§ + ij + i§ + l{ [ig] . This yields 

U t = e AAt U t0 + QT A * e As ds^j F(U t0 ) + £ e A ^BdW s 

+ f e A ^F'(U t0 )(e AAs - I)U tQ ds + 0((At)( 1 + min ( 2 ^)) 
Jt 

for the Taylor expansion corresponding to the S-wood W2. 
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5.1.4. Taylor expansion of order 1 + 2min(7, 6). In examples, we often- 
times have 7 = \ — e, S = \ for an arbitrarily small e G (0, 4) (space-time 
white noise, see Section 5.4) or 7 = ^ — e, <5 = 5 for an arbitrarily small 
e G (0, ^) (trace-class noise). In these cases, the sixth stochastic tree in 
W2 turns out to be the active tree of the lowest order. Therefore, we con- 
sider the stochastic wood W3 = £( 6>1 )W2 (see Figure 5) here. It has order 
1 + 2min(7,<5) (see Section 4.4) and we have 

*(W 3 ) = I° + i? + I° 2 + I{ [ig] + /J [I 2 °] , 

which implies 

U t = e AAt U t0 + ^ e As ds^j F(U tQ ) + jf e^-^Bd^ 

+ f e A ^F'(U to )(e AAs -I)U to ds 
J t 

+ I e A ^F'{U to ) f\ A ^BdW r ds + 0((At) 1+2min ^). 

J to J to 

This example corresponds to the Taylor expansion introduced in the begin- 
ning in Section 3.3. 

5.1.5. Taylor expansion of order 1 + min(37,7 + 5, 26). Since we often 
have 7 < p and 7 < \ in the examples below, the seventh stochastic tree in 
W3 has the lowest order in these examples. Therefore, we consider the Taylor 
approximation corresponding to the S-wood W4 = E(jnW3 (see Figure 6), 
which is given by 

U t = e AAt U t0 + (^j\ As ds^F{U t0 ) + j\ A ^BdW s 

+ f 'e A ^F'(U t0 )(e AAs - I)U t0 ds 
ho 

+ ~ f e A ^ s) F"{U t0 ){^ AAs ~ I)Ut ,(e AAs ~ I)U to )ds 
1 J to 

+ f e A ^F'{U to ) [ S e Ai - s -^BdW r ds 
ho ho 

+ 0((At)( 1+min ^ 37 ' 7+5 ' 2<5 ^). 
It is of order 1 + min(37, 7 + 5, 25), which can be seen in Section 4.4. 
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5.1.6. Taylor expansion of order 1 + 3min(7, 5, |). In the case 7 < ^ 
the 9th, 10th and 12th stochastic tree in W4 all have lower or equal orders 
than the fifth stochastic tree in W4. Therefore, we consider the S-wood W5 = 
-E( 12i i)-E'( 1 o i i)-E(9 i i)W4 (see Figure 7) with the Taylor approximation 

U t = e AAt U t0 + (^ A< e As ds^j F(U t0 ) + j' e A ^BdW s 

+ f 'e A ^F'(U t0 )(e AAs - I)U tQ ds 
Jt 

+ I f e A(t-s) F „ {Uto){{e AAs _ I)Ukv{e AAs _ I)UtQ)ds 

1 Jto 

+ f e A ^F"(U t0 )((e AAs -I)U t0 , f 'e A ^ r) BdW r ) ds 
Jto V Jto J 

+ - f e A ^F"(U to )( f S e A ^BdW r , f e A ^BdW r ^\ ds 

2 Jto \Jto Jt J 

+ t e A{t - s) F'(U tQ ) /% A ( s - r )SdVF" r ds + 0((At)( 1+3min ( 7 ' 5 ' 1 / 3 ))). 

J to J to 



Remark 1. Not all Taylor expansions for general finite-dimensional 
SODEs in [29] are used in practice due to cost and difficulty of computing the 
higher iterated integrals in the expansions. For SODEs with additive noise, 
however, the Wagner-Platen scheme is often used since the iterated integrals 
appearing in it are linear functionals of the Brownian motion process, thus 
Gaussian distributed and hence easy to simulate. A similar situation holds 
for the above Taylor expansions of SPDEs. In particular, the conditional 
distribution [with respect to F'(Ut )] of the expression 

f e A{t ~ s) F'{U t0 ) f e A(s ~ r) BdW r ds 
Jto Jto 

for t £ [to,T] in Section 5.1.4 is Gaussian distributed and, in principle, easy 
to simulate (see also Section 6). 

5.2. Taylor expansions for finite-dimensional SODEs. Of course, the ab- 
stract setting for stochastic partial differential equations of evolutionary type 
in Section 2 in particular covers the case of finite-dimensional SODEs with 
additive noise. The main purpose of the Taylor expansions in this article is 
to overcome the need of an ltd formula in the infinite-dimensional setting. 
In contrast, in the finite-dimensional case, Ito's formula is available and the 
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whole machinery developed here is not needed. Nevertheless, we apply in 
this subsection the Taylor expansions introduced above to stochastic ordi- 
nary differential equations with additive noise to compare them with the 
well-known stochastic Taylor expansions for SODEs in the monograph [29]. 
These considerations are not so relevant in the view of applications, since 
the finite-dimensional case is well studied in the literature (see, e.g., [35] and 
the above named monograph), but more for a theoretical understanding of 
the new Taylor expansions introduced here. More precisely, only in this sub- 
section let H = R rf with d S N be the d-dimensional M-Hilbert space of real 
(i-tuples with the scalar product 

(v, w) = vx ■ w\ H \-v d ■ w d 

for every v = (yi, . . . , Vd) G H and every w = (w±, . . . , Wd) G H. Let also U = 
M m with meN and suppose that (Wt)t£\o,T] is a standard m-dimensional 
Brownian motion. Furthermore, we suppose that the eigenfunctions and the 
eigenvalues of the linear operator —A in Assumption 1 are given by e% = 
(1,0,..., 0) eH, e 2 = (0,1,0,..., 0) £H,...,e m = (0, . . . , 0,1) G H and \ x = 
A2 = • • • = Xd = with X = {1, . . . , d}. So, in this case A is of course a boring 
bounded linear operator with D(A) = H = M. d and Av = for every v £ 
D(A). Furthermore, note that D((k — A) r ) = H = M. d for every r £ R with 
an arbitrary k > 0. The bounded linear mapping i?:R m — > M. d is then a 
d x m-matrix. Due to Assumption 2, the drift term F:R rf — >~R. d is then 
a smooth function with globally bounded derivatives as it is assumed in 
[29]. The initial value xq : Q — > M. d is then simply a J r o/S(R rf )-measurable 
mapping, which satisfies E|xo| p < 00 for every p £ [l,oo). So, the SPDE (1) 
is in that case in fact a SODE and is given by 

(23) dU t = F(U t )dt + BdWt, U = u , 

for t S [0, T] . Now, we apply the abstract Taylor expansions introduced above 
to that simple example. Therefore, note that the parameters in Assumption 
3 are given by 7 = 1 — e and 5 = \ for every arbitrarily small e G (0, 1). First 
of all, we have 

U t = U t0 +B-(W t -W t0 ) + O(At) 

(see Section 5.1.1). Thus, this Taylor approximation corresponds in the case 
of finite-dimensional SODEs to the Taylor approximation for SODEs with 
the multi-index set 

A = {v,(l),(2),...,(m)} 

in Theorem 5.5.1 in [29]. Here, we only mention the multi-index set, which 
uniquely determines the stochastic Taylor expansion in [29] and refer to the 
above named monograph for a detailed description of the stochastic Taylor 
expansions for SODEs there. 
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The exponential Euler approximation in Section 5.1.2 yields 

(24) U t = U t0 + F(U t0 ) -(t-t ) + B- (W t - W t0 ) + 0((At) 3 / 2 ). 

This is nothing else than the corresponding one-step approximation of the 
classical Euler-Maruyama scheme (see Section 10.2 in [29]) and is in the 
setting of [29] given by the multi-index set 

4 = {u,(0),(l),(2),...,(m)} 

in Theorem 5.5.1 there. In that sense, the name of the exponential Euler 
scheme is indeed justified. While in this article, the Taylor approximation 
(24) is obtained via an expansion of the Ij-operators (see Lemma 1 and 
Section 3.1), in [29] the stochastic Taylor approximation (24) is achieved by 
applying Ito's formula to the integrand F(Ut) in the SODE (23). Finally, 
the Taylor approximation in Section 5.1.4 reduces to 



U t 



to I 



U to +F(U t0 )-(t-t ) + B-(W t 
+ F'{U t0 )-B-( f f dW r ds)+0((At) 2 ). 

\Jto Jto J 

The approximation above is nothing else than the one-step approximation 
of the stochastic Taylor approximation given by the multi-index set 



A 



( V ,(0),(l),(2),...,(m), 

(1.0) ,(2,0),...,(m,0), 

(1.1) ,(2,l),...,(m,l), 

, (l,m),(2,l),...,(m,m) , 



in Theorem 5.5.1 in [29]. In [29], it is obtained via again applying Ito's 
formula. 

To sum up, although the method for deriving Taylor expansions in this 
article (Ij-operators) is different to the method in [29] (Ito's formula), the 
resulting Taylor approximations above coincide. 



5.3. Simultaneous diagonalizable case. We illustrate Assumption 3 with 
the case where A and B are simultaneous diagonalizable (see, for example, 
Section 5.5.1 in [6]). This assumption is commonly considered in the liter- 
ature for approximations of SPDEs (see, e.g., Section 2 in [33] or see also 
[24, 32, 34] ). Suppose that U = H and that B : H -> H is given by 

Bv = 6j(e«, w)ej Mv£H, 
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where b{ , i El, is a bounded family of real numbers and , is the 

family of eigenfunctions of the operator A (see Assumption 1). Concerning 
Assumption 3, note that 

[ T \( K - A ye As B\ 2 HS ds 
Jo 

= g( K + A,) 2 ^(^ T e- 2A -^) 



,2kT 



Hex 

JIkT 



^ 6 2 (At + A4) (2 7 -l) (1 _ e -2(A l+K )T ) A 



„2/tT / \ 



for a given 7 > 0, so 



f T \( K - A ye As B\ 2 ls <^ 
Jo 



follows from 

(25) ^2^( K + Xi)^-^<oo 

for a given 7 > 0. In this case, we also have 

(26) f^B^ds^Ct 25 

Jo 

for every t £ [0, 1] with 5 := min(7, 1) and a constant C > 0. 

5.4. Space-time white noise. This example will be a special case of the 
previous one. Let H = L 2 ((0, 1),R) be the space of equivalence classes of 
square integrable measurable functions from the interval (0,1) to R with 
the scalar product and the norm 



(u, v) = J u(x)v(x)dx, \u\ = J J \u(x)\ 2 dx, u,v£H. 
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Let U = H and let B = I:H — > H be the identity operator. In addition, 
assume that a: (0, 1) — >■ K is a bounded measurable function and that the 
operator F :H —> H is given by 

F(v)(x) := (Fv)(x) := a(x) ■ v(x), x G (0, 1), 

for all v G H , which clearly satisfies Assumption 2. Also note that 

F'(v)w = F(w) and . . . , Wi ) = 

for all v,w,w\, . . . ,Wi €: H and all i € {2, 3, . . .}. Furthermore, let A = : 
D(A) C -ff — )• H be the Laplace operator with Dirichlet boundary condition, 
that is, 

oo 

Au = -A n (e n ,u)e n , u £ H, 

71=1 

where 

\ n =ir 2 n 2 , e n (x) = V / 2sin(n7rx), i£ (0,1), 

for each n G X := N. Of course, the e n , n G N, form an orthonormal basis of 
(Assumption 1). Additionally, we choose k = 0. 

Let to = and T = 1. In view of (25), Assumption 3 requires 7 = \ — e for 
every arbitrarily small e > 0. However, instead of (26), we obtain here the 
stronger result S = j, since 

|^ s |hs<C^ 

for every s G (0,1] and a constant C > (see Remark 2 in [32]). Finally, 
let Uo £ H be an arbitrary (deterministic) function in H, which satisfies 
Assumption 4. The SPDE (1) is then given by 



dU t (x) 



d 

T—^Utix) + a(x)U t {x) 
ox z 



dt + dW t , U t (0) = U t (l) = 0, 



with Uq(x) = uo(x) for x G (0, 1) and t G [0, 1]. After considering Assump- 
tions 1-4 for this example, we now present the Taylor approximations in this 
case. Here, e G (0, 4) is always an arbitrarily small real number in (0, 7). 

5.4.1. Taylor expansion of order 1. For an approximation of Ut of order 
one for small t > 0, we obtain 

U t = e At u + f* e A(t - s) dW s + O(At) 
Jo 

(see Section 5.1.1). 
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5.4.2. Taylor expansion of order 4 — e. Here, we have 

U t = e At u + A-\e At - I)Fu + f e A{t ~ s) dW s + 0((At)( 5 / 4 " e )) 

J o 

for an approximation of order | — e (see Section 5.1.2). 

5.4.3. Taylor expansion of order |. In the next step, we obtain 

U t = e At u + A~ 1 (e At - I)Fu + f e A ^~ s) dW s 

Jo 

+ (J e A ^F(e As - I) ds^j u + <3((At) 5 / 4 ) 
for an approximation of order | (see Section 5.1.3). 

5.4.4. Taylor expansion of order | — e. Here, we have 

U t = e At u + A- 1 (e At - I)Fu + e A ^F^ e A ^ dW^j ds 

+ Qf e A{t - s) F(e As - I) ds^j u + jf e A( *~ s) dW s + 0((Atf /2 - £) ) 
for an approximation of order § — £ (see Section 5.1.4). 

5.4.5. Taylor expansion of order X — e. Since F is linear here with F'{y) = 
F and F"(v) = for all v G H , the approximation above is even more of order 
| — e, that is, 

U t = e At u Q + A-\e At - I)Fu + jf e A ^F^ e A ^ dW r ^j ds 

+ Qf e A (*- s )F(e As - I) ds) n + jf e A (*" s ) dW s + 0((Atf^) 
(see Section 5.1.6). 

5.4.6. Taylor of order 2 — e. We also consider the Taylor expansion given 
by the stochastic wood 

^(24,1)^(22, l)^!,!)^^,!)^^, 1)^(16,1)^(15, 1)^(13,1) w 5) 

where W5 is presented in Figure 7. Since F is linear here, we see that the 
corresponding Taylor approximation is the same as in the both examples 
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above, so we obtain 



U t = e At u + A-\e At - I)Fu + jf* e A (*-^ F '(J* e A ^ dW^j ds 

+ (J e A{t - s) F(e As - I) ds^j n + J e A(i ~ s) dW s + 0({Atf- £ 1). 

By further expansions, one can show that this approximation is in fact of 
order 2. 

5.5. Trace class noise. In this subsection, we compute the smoothness 
parameters 7 and 5 in Assumption 3 for the case of trace class noise (see, e.g., 
Sections 4.1 and 5.4.1 in [6]). This assumption is also commonly considered in 
the literature for approximations of SPDEs (see, e.g., [15, 33]). Precisely, we 
suppose that B : U — )■ H is a Hilbert-Schmidt operator, that is, |.B|hs < °°- 
Hence, we obtain 



e As B\ 2 ds < I *(| e As | 2 • \B\la) ds < e 2K \B\ 2 *t 



/ |' -'->lns"' s / II' l ' "Ins. 1 "'" ( "Ins 
Jo Jo 



and therefore 



'J \e As B\ 2 ls ds<e K \B\ US Vi 



for all t 6 [0, 1]. Moreover, we have 

i-T 



I 2 ds 

Ihs as 



[ \(K-A) r e As B\^ s ds< [ \(k - A) r e As \ 2 \B 
Jo Jo 

£\(k - A) r e {A - K)s e KS \ 2 ds^j \B 
£\( K -AY^ A - K >\ 2 d s y 2KT \B\ 2 



2 

HS 



< 



< ( / s~ 2r ds)e 2KT \B\ 2 m <oo 



for all r 6 [0, i). Hence, we obtain 7 = ^ — 6 and S = h for every arbitrarily 
small e € (0, 5) in this situation. Now, we present the Taylor expansions 
from Section 5.1 again in this special situation. 

5.5.1. Taylor expansion of order 1. Here, we have 

U t = e AAt U to + [ e A ^ B dW s + O(At) 
Jt 

for a Taylor approximation of order 1 (see Section 5.1.1). 
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5.5.2. Taylor expansion of order | — e. For a Taylor approximation of 
order | — e (see Section 5.1.2), we obtain 

U t = e AAt U t0 + ^ A ^ As d S ^([/ t J + ^ < e A ( ^ - s ) J B^ + 0((A^)( 3 / 2 - £ )). 

Here and below, e € (0, ^) is an arbitrarily small real number in (0, i). 

5.5.3. Taylor expansion of order | . The Taylor approximation in Section 
5.1.3 reduces to 

U t = e AAt U t0 + {^f\ As ds^F{U t0 ) + f\ A ^BdW s 

+ f* e A ^F'(U t0 )(e AAs - I)U t0 ds + 0((Atf 2 ). 
J t 

5.5.4. Taylor expansion of order 2 — e. Here, we obtain 

U t = e AAt U t0 + e As ds^F(U to ) + f\ A ^BdW s 

+ f e A ^F'{U tQ ){e AAs -I)U tQ ds 

Jto 

+ f e A{t - s) F'{U tQ ) ^ e A{s - r) BdW r ds + 0({Atf-^) 

Jto Jto 

for a Taylor expansion of order 2 — e. This example corresponds to the Taylor 
expansion introduced in Section 5.1.4. 

5.6. A special example of trace class noise. Let H = U = L 2 ((0, l) 3 , R) 
be the space of equivalence classes of square integrable measurable functions 
from (0, l) 3 to R and consider two distinct Hilbert bases ej, i Sl:=N 3 , and 
fi, i €l, in H given by 

ei(x±,X2,x 3 ) = 2 3//2 sin(ii7rxi) sin(i 2 7rx 2 ) sin(i 3 7ra;3) 

and 

fi(x 1 ,x 2 ,x 3 ) 

= C( il _i)C( i2 _ 1 )C(i 3 _ 1 ) COs((il - l)7TXl) COs((i 2 ~ l)ltX 2 ) C0S((« 3 - 1)7TX 3 ) 

for every i = i 2 , is) G T = N 3 and every x = (x\, x 2 , x 3 ) £ (0, l) 3 , where 
c n := v2 for every n G N and cq = 1. Then, consider the Hilbert-Schmidt 
operator B:U ^ H given by 
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for all u 6 U = H . Moreover, let Aj, i G N 3 , be a family of real numbers 
given fry Aj = 7r 2 (z_| + i\ + «§) for all % = (11,12,^3) G N 3 . Finally, consider 
A=(J^ + J| r + J^):.D(A)c.ff->-.ff (with Dirichlet boundary conditions) 

given by 

Av = -\i(ei,v)ei 

igN 3 

for all v G -D(^4), where -D(A) is given by 



^(^ + i| + ii)|(ei,v)p 



Then, the SPDE (1) reduces to 

Q 2 Q 2 Q2 



dUAx) 



dxj + M + oJ ]Ut{x) + Fmx)) 



dt + ^QdW t (x) 



b 2 ^3- 

with f7|d(o,i)3 = for x G (0, l) 3 and t G [0,T]. Assumptions 1-4 are fulfilled 
with S = 5 and 7 = ^—5 for every arbitrarily small e G (0, 5) (see Section 
5.5). The Taylor approximations in that situation are presented in Section 
5.5. 

6. Numerical schemes based on the Taylor expansions. In this section, 
some numerical schemes based on the Taylor expansions in this article are 
presented. We refer to [21, 22, 24, 25] for estimations of the convergence 
orders of these schemes and also for numerical simulations for these schemes. 

For numerical approximations of SPDEs, one has to discretize both the 
time interval [0, T] and the R-Hilbert space H. For the discretization of the 
space H, we use a spectral Galerkin approximation based on the eigenfunc- 
tions of the linear operator A : D(A) C H — ^ H. More precisely, let (Xjv)jveN 
be a sequence of increasing finite nonempty subsets of X, that is, 7^ In C 
T M C X for all N,M G N with N < M and let H N : = span(ei,z G X N ) be 
the finite-dimensional span of |X7v|-eigenfunctions for N G N. The bounded 
linear mappings Pjy ■ H — > Hjy are then given by Pn(v) = ^2 ie .x N ( e i> v ) e i f° r 
every v G H and every N G N. 

6.1. The exponential Euler scheme. Based on the Taylor approxima- 
tion in Sections 3.2 and 5.1.2, we consider the family of random variables 
Y*' M : n -»• H N , k = 0, 1, . . . , M, N, M G N, given by Y Q N ' M = P N (u Q ) and 



Y, 



N,M 
k+1 



e ATjMyN,M 



+ 



(27) 



T/M 







e As ds )(P N F)(Y< 



N,M\ 



+ Pn 



(fc+l)T/M 



kT/M 



e A((k+l)T/M-s) BdW 
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for every k = 0, 1, . . . , M — 1 and every N,M £~N. This scheme is introduced 
and analyzed in [24]. As already mentioned, it is called the exponential Euler 
scheme there. 

In the setting of deterministic PDEs, that is, in the case B = 0, this 
scheme reduces to 

Y N,M = e AT/M Y N,M + f^/M ^ X ^^N^ 

for every k = Q,l,...,M — 1 and every JV,MgN. This scheme and similar 
schemes, usually referred as exponential integrators, have for deterministic 
PDEs been intensively studied in the literature (see, e.g., [3-5, 16-19, 27, 
30, 31, 36]). Such schemes are easier to simulate than may seem on the first 
sight (see [4]). In the stochastic setting, we refer to Sections 3 and 4 in [24] 
for a detailed description for the simulation of the scheme (27), in particular, 
for the generation of the random variables used there. 



6.2. The Taylor scheme indicated by Section 5.1.3. In view of Section 



5.1.3, we obtain the Taylor scheme Y. ' : U -> H N , k = 0, 1, . . . , M, N, M £ 



N, given by Y^' M = Pn(uq) and 



Y N,M = e AT/M Y N,M + ^ & As ^ ( y N, 

r(k+l)T/M 

+ / e M ^ T l M - s \P N F')(Y h 

J kT/M 



k 

I kT/M 

x{{e A(s-kT/M)_ I)Y N,M )ds 
/ r {k+l)T/M \ 

+ P N ( e A ^ T l M -^BdWs) 

\JkT/M J 

for every k = 0, 1, . . . , M — 1 and every N, M G N. 

6.3. The Taylor scheme indicated by Section 5.1.4- The Taylor approx- 
imation in Section 5.1.4 yields the Taylor scheme Y, ,M :Q — > Hn, k = 
0, 1, . . . , M, N, M e N, given by Y Q N ' M = P N (u ) and 



Y N,M = e AT/M Y N,M + f£ /M e As A 



k^) 



(k+l)T/M 

+ / e A{{k + l)T/M~s) { p NFl){Y N,M , 

kT/M 

x{{e Ms-kT/M)_ I)Y N,M )ds 
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r (k+i)T/M 

+ / e M{k + l)T/M-s) { p NF , ){ yNM ) 
JkT/M 

x (p N e A( - s - u) B dW^j J ds 

/ r (k+l)T/M \ 

+ P N ( e A ^ T l M -^BdW s ) 

\JkT/M J 



' kT/M 

for every k = 0, 1, . . . , M - 1 and every iV, M G N. 

6.4. ^4 Runge-Kutta scheme for SPDEs. In principle, we can proceed 
with the next Taylor approximations and obtain numerical schemes of higher 
order. These schemes would however be of limited practical use due to cost 
and difficulty of computing the higher iterated integrals as well as the higher 
order derivatives in the Taylor approximations. Therefore, we follow a dif- 
ferent approach and derive a derivative free numerical scheme with simple 
integrals — a so called Runge-Kutta scheme for SPDEs. We would like to 
mention that this way is the usual procedure for numerical schemes for differ- 
ential equations: Taylor expansions and their corresponding Taylor schemes 
provide the underlying theory for deriving numerical schemes, but are rarely 
implemented in practice. Instead of these Taylor schemes other numerical 
schemes, which are easier to compute but still depend on the Taylor expan- 
sions such as Runge-Kutta schemes or multi-step schemes (see, e.g., [8] for 
details) are used. 

To derive a Runge-Kutta scheme for SPDEs, we consider the Taylor ap- 
proximation in Section 5.1.4 (see also the Taylor scheme above) from to 

^ k+ M T ana - obtain 

U {k+1)h ^e Ah U kh +( [ h e As ds]F(U kh )+ I * + * " e WM B ^ 

\J0 / Jkh 

+ [ k+lh e A ^ k+1 ^F'(U kh )((e A ^- kh ^-I)U kh )ds 

Jkh 

r(k+l)h 



+ / e A ^ k+1 ^F'{U kh )( e A ^BdW r )ds 

Jkh \Jkh J 



and hence 



r(k+l)h 

U, k+1)h « e Ah U kh + he Ah F(U kh ) + / e A « k+1 ^BdW f 

Jkh 
/ 1 Ak+l)h 

+ he Ah F'{U kh ) (- J ( e Ms-kh) _ I)Ukh ds 
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/ 1 p(k+l)h ps \ 

+ he Ah F'(U kh ) - / / e A ^B dW r ds 

\ n Jkh Jkh J 

e Ah U kh + he Ah F{U kh ) + e A ^ k+1 ^ B dW s 

Jkh 



+ he Ah F'(U kh ) 



-i r (k+l)h r kh 

h 



kh 

/ i p(k+l)h pa 

+ he Ah F'{U kh )(- / e^-^BdWrds 

\h Jkh Jkh j 

with h := -21 for k = 0, 1, . . . , M - 1 and M G N. This yields 

~ e^Cfa, + he Ah F(U kh + Zjf) + e A ^ k+1 ^BdW s 

Jkh 

with the random variables 

i p(k+l)h pkh 
Z^ = =- ( e A{s-kh)_^ / e A{kh-r) BdWrds 

h Jkh Jo 

1 p(k+l)h ps 

+ - I I e A( - s -^BdW r ds 



pKIl 

e A(s-kh)_ I j / e A(fch-r) Bdwi-da 
./0 



fc/i ■/ fc/i 

for k = 0,1, . . . ,M — 1 and M 6 N. The corresponding numerical scheme 
Y k NM H N ,k = 0,1,..., M,N,MeN, is then given by Y N ' M = P N {u ) 
and 

yjJMf = e AT/M Y N,M + ^AT/M^^M + p^M^ 

/ p(k+l)T/M \ 
\JkT/M J 

for every = 0, 1, . . . , M — 1 and every iV, M G N. This Runge-Kutta scheme 
for SPDEs is introduced and analyzed in [21]. Under non-global Lipschitz 
coefficients of the SPDE, it is analyzed in [22]. Note that the random vari- 
ables occurring in the scheme above are Gaussian distributed and therefore 
easy to simulate (see also Remark 1 and [21, 22] for details). More pre- 
cisely, in the case of one-dimensional stochastic reaction diffusion equations 
with space-time white noise it is shown in the articles cited above that this 
scheme converges with the overall order \ — with respect to the number of 
independent standard normal distributed random variables and the number 
of arithmetical operations used to compute the scheme instead of the overall 
order i of classical numerical schemes (see, for instance, [7, 12, 13, 42]) such 
as the linear implicit Euler scheme. 
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7. Proofs. 

7.1. Proo/s of (6) and (7). 

Lemma 1. Let Assumptions 1~4 be fulfilled and let i G N be given. Then, 
we have 

i 1 ,=/ 1 + / 1 M/ ] + / 1 M/ 1 ,] + / 1 M/ 2 °] 

and 

i\, [ 91 ,. ..,#]= . ..,#]+ if; +1) [/ ,5i, • ■ .,<*] 

+/{: +1) [/ 1 ,, 5l ,...,< 7i ]+/it +1) [/ 2 °, 5l ,..., ffi ] 

/or a// g 1 ,...,g i G P. 

Proof. We begin with the first equation. Since we have 

F(U S ) = F(U t0 ) + f F'(U t0 + r(U s - U to ))(U s - U to )dr 
Jo 

= F(U to )+ f F'(U to +rAU s )(AU s )dr 
Jo 

= F(U t0 )+ f F'(U t0 +rAU 8 )(I$(s))dr 
Jo 

+ / F'(U t0 + rAU s )(I^(s))dr+ C F'(U to +rAU s )(I$(s))dr 
Jo Jo 

for every s G [to,T] due to the fundamental theorem of calculus and (5), we 
obtain 



ij.(t)= / e A ^F(U s )ds 

J t 

= f e Mt- s ) F{Uto)ds+ f\A{t-s)( f 1 F >(u tQ+ rAU s ){I*(s))dr\ds 
Jt Jt \Jo J 

+ £ e A ^ ( jf F'(U t0 + rAU.Xq. (a)) dr^j ds 

+ ft^^ (i F '( U to+r&U s )(I$(s))dr^ ds, 
which implies 

Jf. (t) = /? (t) + I\» [I° ] (t) + I}. [ij. ] (t) + /J. [I 2 °] (0 
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for all t G [io>2~1- Moreover, we have 

F®(U t0 + rAC/ s )( ffl ( S ), . . ., 9i (s)) { \ dr 
o ^ _i J ! 
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-F®{U to +rAU s ){ 9l (s),..., 9i ( S )) 



(l-rf 



r=l 



r=0 



+ / F^ i+1 \U to +rAU s )(AU s , 9l (s),..., 9i (s)) 



(1-rf 



dr 



l -F^{U tQ )( gi (,),..., 9i (s)) 



+ / F( i+1 \u t0 + rAU s )(AU s , 9l (s), . . ., 9i (s))^A- dr 
Jo 



for all s G [to; ^1 and all gi, . . . , 9 i G V due to integration by parts and there- 
fore, we also obtain the second equation. □ 

7.2. Proof of Theorem 1. For the proof of Theorem 1, we need the fol- 
lowing lemma. 

Lemma 2. Let X : [to, T]xO-> [0, oo) be a predictable stochastic process. 
Then, we obtain 



X s ds 



t 



< / \X s \ L rds 



for every t G [toj^] an d every r G [1, oo), where both sides could be infinite. 



Proof. First, we consider the case, where Xt < C is bounded by a 
constant C > for all t G [0,T]. Here, we have 



E 



X s ds 



to 



E 



to 



X,, du 



to 



(r-1) 



X, 



< / A'.. ,.. ds[ E 



to 



((r-l)/r) 



for every i G [io,^] and every r G [l,oo) due to Holder's inequality. Since 

X u du 



E 



to 



< OO 



is finite for every t G [to,T] and every r G [l,oo) due to the boundedness 
of X : [to,T] x fi — > [0,oo), we obtain the assertion. In the general case, we 
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can approximation the stochastic process (^Q)tg[o,Tl by bounded processes 
(X t N ) te[0 , T] for N G N given by 



X? :=mm(N,X-, 



for all t G [0,T] and all Af£N. This shows the assertion. □ 

We also need the Burkholder-Davis-Gundy inequality in infinite dimen- 
sions (see Lemma 7.7 in [6]). 



Lemma 3. Let X : [to, T] x — > HS(U, H) be a predictable stochastic pro- 
is with E f fo \X S \^ S < 00 • Then, we obtain 

1/2 



t 

X s dW s 

to 



<p[ I \\X s \us\lp ds 



LP \J t 

for every t G [io^] an d every p G [l,oo), where both sides could be infinite. 

In view of the definitions of the mappings <I> and Theorem 1 immedi- 
ately follows from the next lemma. For this, the subset ST' C ST of stochas- 
tic trees given by 

ST' := jt = (t',t") G ST|VA: G nd(t) : ((3/ G nd(t):t'(Z) = k) 

=>t?(k)e{l,l*}))} 

is used. 

Lemma 4. Let t = (t',t") G ST' be an arbitrary stochastic tree in ST'. 
Then, for each p > 1, there exists a constant C p > such that 

(E[|$(t)(t)H) 1/p <C p -(t-to) ord(t) 

holds for all t G [to,T], where C p is independent oft and to but depends on 
p,t,T and the SPDE (1). 

Proof. Due to Jensen's inequality, we can assume without loss of gen- 
erality that p G [2, oo) holds. We will prove now the assertion by induction 
with respect to the number of nodes l(t) G N. 

In the base case Z(t) = 1, we have $(t) = L®,,,^ by definition. Hence, we 
obtain 

MW)\lp = \4, {1) (t)\ LP = \L° (t)\ LP = \(e AAt - L)U t0 \ LP 
< |(« - A)^(e AAt -L)\-\(k- Ayu t0 \ LP 



TAYLOR EXPANSIONS FOR SPDE'S 



35 



|(k - A)~~< (e {A - K)M - e~ KAt )\ ■ e nAt 



< \ (K-Ay~<( e ( A - K)At 



e - KAt )\-e KT - 



\( K -Ayu t0 \ LP 
sup \{K-Ayu s \ LP 

0<s<T 



<C p -\(K-A)'^(e {A - K)At 



-nAt\ 



and therefore 

MW)\lp < C p • (|(« - A)^(e^- K ) A * - 7)| + |(« - A)-^(J - e- KA *)|) 

< C p • (|(« - A)^{e^ A -^ At - I)\ + |(I - e- KA ')|) 

< C p ■ (At) 7 + C p ■ (At) < C p ■ (At) 7 

for every t E [to,T] in the case t"(l) = 0. Here and below, C p > is a con- 
stant, which changes from line to line but is independent of t and to- More- 
over, by Lemma 2, we obtain 

r-t 



MW)\ 



LP 



LP 



to 



< 



< 



Iv'(i)(t)\ LP = \lUt)\ 

(\e A ^F(U s )\ LP )ds<C p 



^- s )F{U s )ds 



LP 



to 



\F(U a )\LP ds 



to 



Cp-^j\l + \U s \ LP )ds^<C p -(At) 



for every t E [to,T] in the case t"(l) = 1* and 

\m(t)\LP = \i° t , [1) m L p = \ii(t)\tp 

" e A(ts) F( U to )ds 



to 



LP 



< 



< 



(\e A ^F(U t0 )\ LP )ds<C p 



to 



\F(U t0 )\ LP ds 



to 



Cp-^j\l + \U t0 \ LP )ds^<Cp-(At) 



for every t E [to, T] in the case t"(l) = 1. Finally, due to Lemma 3, we obtain 



Mm)\LP = \i t , {1) (t)\ LP = \i° 2 (t)\ 



LP 



<c p - 



<c n 



le A(t-s) Bl ,2 



to 



HSI LP d s j 



e A (ts) BdWs 

\ 1/2 



to 



LP 



(At) 



1/2 
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\S 



< C p ■ (At) 6 

for every t G [t ,T] in the case t"(l) = 2. This shows \<S>(t)(t)\ LP < C p ■ 
(Ai) ord W for every t G [t ,T] in the base case l(t) = 1. 

Suppose now that l(t) G {2,3, . . .}. Since t = (f ,t") € ST', we must have 
t"(l) G {1,1*}. Let ti,...,t„ G ST' with n G N be the subtrees of t. Note 
that ti, . . . ,t n are indeed in ST'. Then, by definition, we have 

$(t)(t)=i?, (1) [$(ti),...,*(t n )](t) 

for every t G [ioj^l- Therefore, by Lemma 2, we obtain 

|<D(t)(t)| LP = |/ t "»(i) [*(ti), • • • , *(t„)](t)| = |/f. [*(ti), . . • , *(t n )](i)| LP 



f'e^t-^f f 1 F (n) {Uto+rAUs) 
t \Jo 



x (*(ti)(a), . . . , *(t n )(s))^ >— dr ) ds 



(n-1)! 



< 



C P - T ! F^ n \u t0 +rAU a 
Jt Jo 



LP 



x(gKti)(s),...,5Ktn)(s)) (1 ^"/ dr 



(n-1)! 



LP 




to JO 



x(*(ti)(a),...,*(t„)(s)) 



(1-r) 



n-1 



(n-1)! 



LP 



< 



< 



C P - / / \F( n \U t0 +rAU s )($(t 1 )(s),...,$(t n )(s))\ LP drds 

J to JO 

C P - f f 1 \mt 1 )(s)\---Mt n )(s)\\ LP drds 

J t JO 



and hence 



Mm)\ LP <C p - / ||*(ti)( S )|.-.|*(t„)(a)||ipda 
■/ to 

< C v(yf t |*(ti)(s)|iP« • • • |$(t n )(s)| LP « cfc 
< C p • (^j\(As) ord{tl) ■ ■ ■ (As) ord{tn) )ds 



< C p • (At)( 1 + ord ( t i)+- +ord ( t ™)) = C p ■ (At) 



ord(t) 
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for every t E [to,T] in the case t"(l) = 1*, since Z(ti) , . . . ,l(t n ) < Z(t) — 1 and 
we can apply the induction hypothesis to the subtrees. A similar calculation 
shows the result when t"(l) = 1. □ 

7.3. Properties of the stochastic convolution. 

Lemma 5. Let Assumptions 1 and 3 be fulfilled. Then, there exists an 
adapted stochastic process O : Q — > C([0, T],H) with continuous sample paths, 
which is a modification of the stochastic convolution J e A ^~ s ^ B dW s , ZE 
[0,T], that is, we have 



for all t E [0,T]. 

Proof. Let Z : [0, T] x — > H be an arbitrary adapted stochastic pro- 
cess with 



Jo 

for all t E [0,T]. Due to Assumption 3, such an adapted stochastic process 
exists. Moreover, Z : [0, T] x 17 — > H is centered and square integrable with 



for all t E [0, T] . Let 8 := min(<5, 7), p > 2 and < t\ < t<i < T be given. Then, 





P-a.s. 




we have 





P-a.s. 



and 
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due to Lemma 3. Hence, due to Assumption 3, we obtain 
\Zt 2 — Z tl \lp <C(t2 — h) 5 

+ C\(k- A)^(e A ^ - I)\ QT 1 1 (« - Aye As B\ 2 ns ds^j ^ 

/ f T \l/2 

< C(t 2 - h) s + C(t 2 - hVU |(« - A)^ e As i?|^ s d s J 

< C(t 2 - ti)" 5 + C(t 2 - ti) 7 < C(* a - h) e , 

where C > is a constant changing from line to line. Since > is greater 
than zero and since p > 2 was arbitrary, there exists a version of (Zt)t^\o,T] 
with continuous sample paths due to Kolmogorov's theorem (see, e.g., Chap- 
ter 3 in [6]). □ 

Acknowledgments. We strongly thank the anonymous referee for his 
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